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Abstract. Local analysis of the two dimensional Navier-Stokes equations is 
used to obtain estimates on the energy and enstrophy fluxes involving Taylor 
and Kraichnan length scales and the size of the domain. In the framework 
of zero driving force and non-increasing global energy, these bounds produce 
sufficient conditions for existence of the direct enstrophy and inverse energy 
cascades. Several manifestations of locality of the fluxes under these conditions 
are obtained. All the scales involved are actual physical scales in M.'^ and no 
homogeneity assumptions are made. 



1. INTRODUCTION 

Following the groundbreaking ideas of Kolmogorov [15l[T6l[T4], Batchelor, Kraich- 
nan and Leith [2j [3l [171 UHl ED] established the foundations of empirical theory of 
2D turbulence (BKL theory). One of the main features of the BKL theory is the 
existence of enstrophy cascade over a wide range of length scales, called the inertial 
range, where the dissipation effects are dominated by the transport of enstrophy 
from higher to lower scales. In contrast to the 3D turbulence, the energy in 2D 
case is cascading toward the larger scales, a phenomenon referred to as the inverse 
energy cascade. Direct enstrophy and inverse energy cascades have been observed 
in physical experiments (albeit certain difficulties exist in generating a purely 2D 
turbulent flow), but theoretical justification of these phenomena using equations 
of fluid motion, and in particular, the Navier-Stokes equations (NSE), remains far 
from being settled. Technical complexity of the NSE makes it difficult to estab- 
lish the conditions under which such cascades can occur. In the 2D case, the NSE 
possess a number of useful regularity properties (unlike the 3D case for which the 
global regularity is an open problem). However, the dynamical complexity of the 
NSE makes a detailed study of their long time behavior a diflicult enterprise. Under 
certain conditions, existence of the global attractors of high fractal and Hausdorff 
dimensions has been established for the 2D NSE; moreover, it is believed that these 
attractors become chaotic (although the proof is elusive) . For an overview of various 
mathematical models of turbulence and the theory of the NSE, see, e.g., [T2l[T3l[T0] 
and [211 mEH, respectively. 

Most rigorous studies of 2D NSE turbulence have been made in Fourier settings. 
In particular, in |llj the framework of space-periodic solutions and inflnite-time 
averages was used to study main aspects of the BKL theory, including establishing a 
sufficient condition for the enstrophy cascade. This condition, involving Kraichnan 
length scale, is akin to our condition (|4.11l) obtained in section 4. In contrast to 
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[llj . our goal was to work in physical space and with finite-time averages, dealing 
with actual length scales in rather than the Fourier wave numbers. 

In this paper we extend to the 2D case the ideas introduced in [8^ to establish the 
existence of the energy cascade and space locality of the flux for the 3D NSE. There, 
one of the difficulties was the possible lack of regularity, which led us to using the 
framework of suitable weak solutions ([211I1])- In 2D, the difficulties lie in the need 
to work with higher-order derivatives in the case of the enstrophy cascade, as well 
as in dealing with a rather complex phenomenology of the 2D turbulence. 

Despite these differences, the basic setting for studying energy and enstrophy 
transfer in physical scales remains the same in both 3D and 2D case. We utilize 
the refined cut-off functions to localize the relevant physical quantities in physical 
space and then employ ensemble averages satisfying certain optimality conditions 
together with dynamics of NSE to link local quantities to global ones (see [8j for a 
detailed discussion of our physical scales framework). 

We restrict our study to a bounded region, a ball, in R^, and consider the case 
of short-time or decaying turbulence by setting the driving force to zero. Thus, in 
contrast to infinite-time averages used in jllj . we use averages over finite times. 
The time intervals considered here depend on the size of the domain as well as the 
viscosity (see (|4.2p ). The spatial ensemble average is taken by considering optimal 
coverings of the spatial domain with balls at various scales. Also, to exclude the 
situations of the uniform growth of kinetic energy without any movement between 
the scales we restrict our study to physical situations where the kinetic energy on 
the (global) spatial domain is non- increasing, e.g., a bounded domain with no-slip 
boundary conditions, or the whole space with either decay at infinity or periodic 
boundary conditions. 

The paper is structured as follows. In section 2 we provide a brief overview of 
the 2D NSE theory, noting the relevant existence and regularity results. We also 
point out important differences between 2D and 3D NSE, and how these difficulties 
are reflected in the differences between 2D and 3D turbulence. 

Section 3 introduces the physical quantities of energy, enstrophy, and palinstro- 
phy, as well as energy and enstrophy fluxes adopted to our particular settings. We 
also define the ensemble averages to be used throughout the paper. 

The main result of section |4] is a surprisingly simple sufficient condition for the 
enstrophy cascade ()4.11|) . according to which the averaged enstrophy flux toward 
the lower scales is nearly constant over a range of scales. This condition, involving 
the Kraichnan scale and the size of the domain, is reminiscent of the Poincare 
inequality on a domain of the corresponding size (see Remark 14.21) . Moreover, the 
condition in hand would be easy to check in physical experiments as the averages 
involved are very straightforward. 

Section 5 commences a study of inverse energy cascade in physical space. The 
existence of such cascades in the 2D NSE solutions remains an open question. 
Several partial results exist; in particular, in the space-periodic setting the energy 
flux is oriented towards lower (Fourier) scales in the region below the scales of the 
body force ([H]), but existence of the cascade could not be established. In contrast, 
[1] provides a condition for the inverse energy cascades inside spectral gaps of the 
body force. We prove that if the global Taylor scale is dominated by the linear size 
of the domain, then the averaged energy flux is constant over a range of large scales 
and is oriented outwards (see Theorem 15. ip . 
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The second part of the paper concerns locality of the energy and enstrophy fluxes. 
Similarly to the 3D turbulence ([23]), it is believed that the energy and enstrophy 
fluxes inside the inertial ranges of the 2D turbulent flows depends strongly on the 
flow in nearby scales, the dependence on lower and much higher scales being weak. 
The theoretical proof of this conjecture remained elusive. The first quantitative 
results on fluxes were obtained by early 70's (see fT9\). Much later, the authors in 
[22] used the NSE in the Fourier setting to explore locality of scale interactions for 
statistical averages, while the investigation in [5] revealed the locality of filtered en- 
ergy flux under the assumption that the solutions to the vanishing viscosity Euler's 
equations saturate a defining inequality of a suitable Besov space (a weak scaling 
assumption). A more recent work [5i provided a proof of the quasi- locality of the 
energy flux in the Littelwood-Paley setting. 

In section 6 we obtain several manifestations of the locality of both energy and 
enstrophy fluxes in the physical space throughout the inertial ranges. In particular, 
considering dyadic shells at the scales 2''R {k an integer) in the physical space, we 
show that both ultraviolet and infrared locality propagate exponentially in the shell 
number k. 

To the best of our knowledge, the condition ()4.11|) is presently the only condition 
(in any solution setting) implying both the existence of the inertial range and the 
locality of the enstrophy flux. The same is true for the relation (|5.19p which implies 
both inverse energy cascade and energy flux locality in the physical scales of the 
2D NSE. Finally, we point out that our approach is valid for a variety of boundary 
conditions (in particular, the no-slip, periodic, or the whole space with decay at 
infinity); moreover, it does not involve any additional homogeneity assumptions on 
the solutions to the NSE. 

2. PRELIMINARIES 

We consider two dimensional incompressible Navier-Stokes equations (NSE) 

^2 1) x) - z/Au(<, x) + [nit, x) • V)u(i, x) + Vp{t, x) = 

V • u(t,x) ^ , 

where the space variable x = (a:;i, X2) is in and the time variable t is in (0, 00). 
The vector- valued function u = (ui, U2) and the scalar- valued function p represent 
the fluid velocity and the pressure, respectively, while the constant v is the viscosity 
of the fluid. 

Under appropriate boundary conditions this system admits a unique solution (see 
[25) . [6]), which is analytic in both space and time. For convenience, we generally 
assume no-slip boundary conditions on a bounded domain 

(2.2) vl\q^ =0, Vt bounded in 

(although the results hold for the other physical boundary conditions which imply 
smoothness and non-increasing global energy j^^ |up). 

Thus, if e 2?((0,oo) x f2), > 0, where be an open connected set in M^, 
then multiplying NSE by 0u and integrating by parts we obtain the local energy 
equation 
(2.3) 

2i^ / / |V«)Up(/)dxdt= // \vL\'^{dt(j) + vl^(l))dyidt+ // (|up + 2p)u • V(/)dxdt 
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where l?((0,cx)) x ft) denotes the space of mfinitely difFerentiable functions with 
compact support in (0, oo) x fi. 

We also consider the vorticity form of the 2D NSE by taking the curl of (I2.ip 
viewed as a 3D equation with the third component zero, 

d 

(2.4) —oj - i^Auj + (u • V)w = 0, 

where w = V x u (with the convention u — {ui,U2, 0) and uj — (0, 0, w)). 

Note that for the full 3D NSE (12.41) would contain the vortex-stretching term 
(w • V)u. 

Multiplying (|2.4p with 0a; yields the local enstrophy equation, 

(2.5) 2i^yy" |V«)wp(/)dxdi = JJ \u}\'^{dt(l) + vAcj)) dxdt + jj |a;pu • V(/)dxdt . 

We will make the following assumptions on the domain Vt and test functions 4>. 
First, we assume there exists i?o satisfying 

(2.6) i?o > such that B(0, 3i?o) C ^ 

where -8(0, 3i?o) represents the ball in M? centered at the origin and with the radius 
3i?o. 

Next, let l/2<6 < 1. Choose e 2?(B(0, 2i?o)) satisfying 

(2.7) < ^0 < 1, V'o = 1 on B(0, Rq), —g- < — , < -5 ■ 

For a T > (to be chosen later), xo G B{0,Ro) and < i? < i?o, define </) = 
'/'xo,T,_R(i, x) = f]{t)ilj{x.) to be used in (j2.3p and (j2.5p where 77 = ?7T(i) and ijj = 
^xo,-r(x) are refined cut-off hmctions satisfying the following conditions, 



(2.8) 77eI?(0,2T), 0<7/<l, 77 = lon(r/4,5T/4), ^< 

if B{xo,R) C -8(0, i?o), then V € V{B{xq,2R)) with 
(2.9) 

0<^<^o, ^ = loni?(xo,i?)ni3(0,i?o), ^<§- 

and ff B{xq,R) B{0,Rq), then V S V{B{0,2Ro)) with -0 = 1 on B{xq,R) n 

B{0,Rq) satisfying, in addition to (|2.9I) . the following: 

(2.10) 

= 00 on the part of the cone in centered at zero and passing through 
S'(0, -Ro) n -8(xo, R) between S{0, Ra) and s(0, 2i?o) 

and 

= on B{0, i?o) \ -8(xo, 2it!) and outside the part of the cone in 
(2.11) centered at zero and passing through S{0, Rq) n S(xo, 2i?) 
between 5'(0, i?o) and S{0,2Ro). 

Figure [T] illustrates the definition of in the case -8(xo,i?) is not entirely con- 
tained in -8(0, i?o)- 
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!/■ = l/'O 



region containing tlic rest of supp(^/j) 




Figure 1. Regions of supp(V') in the case B{x.o,R) 't B{0,Ro). 



Remark 2.1. The additional conditions on the boundary elements (|2.10p and 
(|2.11l) are necessary to obtain the lower bound on the fluxes in terms of the same 
version of the localized enstrophy E in Theorems 14. Il and l6. 21 (see Remarks 14. 31 and 

EM- 

3. Averaged enstrophy and energy flux 

Let xo G B{0, Rq) and < i? < i?o- We define the localized versions of energy, 
e, enstrophy, E, and palinstrophy, P at time t associated to i?(xo, R) by 

(3.1) exo,flW 



(3.2) E^„Mt) = I 



and 
(3.3) 



1' '2(2(5-1 



idx 



In the classical case, the total - kinetic energy plus pressure - flux through the 
sphere S{xq, R) is given by 

J (i|up +p) u • nds = J ((u- V) u + Vp) •u( 

where n is an outward normal to the sphere ^(xo, R). Similarly, the enstrophy flux 
is given by 

J ^\u}\'^ u ■ n ds — J {u ■ V)uj ■ uj dx . 
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Considering the NSE localized to i?(xo, R) leads to the localized versions of the 
aforementioned fluxes, 



where (p = rjtp with ry and -0 as in (12.8112. 9p . Since -0 can be constructed such that 
Wcj) = rjVil) is oriented along the radial directions of -B(xo, R) towards the center of 
the ball Xq, ^*xo,-R and ^'xo.fl can be viewed as the fluxes into B{xq, R) through the 
layer between the spheres ^(xo, 2R) and ^(xo, R) (in the case of the boundary ele- 
ments satisfying the additional hypotheses (I2.10p and (|2.1ip , ip is almost radial and 
the gradient still points inward). In addition, (|2.3p and (|2.5p imply that positivity 
of these fluxes contributes to the increase of Cxq^r and -Exo.i?,! respectively. 

Note that the total energy flux $xo,i?, consists of both the kinetic and the pressure 
parts. Without imposing any specific boundary conditions on it is possible that 
the increase of the kinetic energy around Xq is due solely to the pressure part, 
without any transfer of the kinetic energy from larger scales into -B(xo, R) (see [8 ). 
As we mentioned in the introduction, under physical boundary conditions, like (j2.2l) . 
the increase of the kinetic energy in _B(xo,-R) (and consequently, the positivity of 
^xo,r) implies local transfer of the kinetic energy from larger scales simply because 
the local kinetic energy is increasing while the global kinetic energy is non-increasing 
resulting in decrease of the kinetic energy in the complement. This is also consistent 
with the fact that in the aforementioned scenarios one can project the NSE to 
the subspace of divergence-free functions effectively eliminating the pressure and 
revealing that the local flux $xo,i?, is indeed driven by transport /inertial effects 
rather than the change in the pressure. 

Henceforth, following the discussion in the preceding paragraph, in the setting of 
decaying turbulence (zero driving force, non-increasing global energy), the positivity 
and the negativity of 'I'xo.-R and 5'xo,-R will be interpreted as transfer of (kinetic) 
energy and enstrophy around the point Xq at scale R toward smaller scales and 
transfer of (kinetic) energy around the point Xq at scale R toward larger scales, 
respectively. 

For a quantity Q — 6x,_r(^), t g [0, 2T] and a covering {i?(xi, i?)}i^i_„ of 
B{0,Rq) define a time-space ensemble average 



(3.4) 




and 



(3.5) 




(3.6) 




Denote by 



(3.7) 



sr = (ex,_R.(*))-R. I 



(3.8) 





(3.9) 



Pr - {PxMt))R 



(3.10) $fl = (*x,fl(t))fl , 
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and 

(3.11) - (^'x,fl(i))fl , 

the averaged localized energy, enstrophy, palinstrophy, and inward-directed energy 
and enstrophy fluxes over balls of radius R covering i?(0, Rq)- 

Also, introduce the time-space average of the localized energy, enstrophy and 
palinstrophy on B{0,Ro), 

(3.12) eo = M l^eoMoit)dt=^^ 11 ^\u\m'-U^dt 



(3.13) 



and 
(3.14) 

where 

(3.15) Mt,^) = vit)M^) 

with tpo defined in (|2.7p . 

Finally, define Taylor and Kraichnan length scales associated with B{0, Rq) by 

(3.16) 



ha 



and 

(3.17) ao . p 

To obtain optimal estimates on the aforementioned fluxes we will work with 
averages corresponding to optimal coverings of B{0,Rq). 

Let Ki,K2 > 1 be absolute constants (independent of R,Ro, and any of the 
parameters of the NSE). 

Definition 3.1. We say that a covering of B{0,Rq) by n halls of radius R is 
optimal if 

(3.8, (^)%„<.,(«.^' 

(3.19) any x G B{0,Ro) is covered by at most K2 balls B{x.i,2R) . 

Note that optimal coverings exist for any < R < Rq provided Ki and K2 are 
large enough. In fact, the choice of Ki and K2 depends only on the dimension of 
e.g, we can choose Ki = K2 = 8. 

Henceforth, we assume that the averages are taken with respect to optimal 
coverings. 

The key observation about these optimal coverings is contained in the following 
lemma. 
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Lemma 3.1. If the covering {-B(xi, of B{0, Rq) is optimal then the aver- 

ages CR, Er, and Pr satisfy 

-^eo <eR< K2eo , 
(3.20) -^^Eo < Er < K2E0 (^^^E'o<E'j,<K2E'^^ , 

-^Po <Pr< K2P0 ■ 

Proof. Note that since the integrand is non-negative, using p.l9p and the lower 
bound in p.lSp we obtain 

Next, we use the upper bound in (13.181) and the non- negativity of the integrand to 
bound cr from below. 



TR^n^JJ 2 - TR^nJJ 2 ^° 

1 1 1 f f \uW2s-1 _ 1 



arriving at the first relation of p.20p . The other two relations are proved in a 
similar manner. 

□ 



Note that the lemma above shows that for the the non-negative quantities, like 
energy, enstrophy, and palinstrophy, the ensemble averages over the balls of size R, 
S-R; Er, and Pr are comparable to the total space-time average. This is not so for 
the quantities that change signs, like the energy and enstrophy fluxes. In fact ^r 
and "^R provide a meaningful information as to energy and enstrophy transfers into 
balls of size R. Positivity of 'S'r, for example, implies that there are at least some 
regions of size R for which the enstrophy flows inwards. 

Moreover, note that the space-time ensemble averages of energy, enstrophy, and 
palinstrophy that correspond to these optimal coverings (over finite number of balls) 
are equivalent to the uniform space-time average. We define the uniform space-time 
average of = 6x,i?(^) as 

2T 

(3.21) / J ^Q^:Rit)d^dt; 

° B{0.,Ro) 

thus we have the following uniform averages of energy, enstrophy, palinstrophy and 
fluxes in regions of size R: e% EJ^ {E'^), P^, $^ and 
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Lemma 3.2. The following estimates hold 
^eo < < 4^60 , 

(3.22) li?o < El < A'Eo ^E', < E'^ < 4^E', ) 



.2^0 <Pr< ^^Po 



1 

22 

Proof. We will prove the first relation in p.22[) . the others follow in a similar way. 
Note that the definition of uniform average applied to the energy ex,_R(i) yeilds 

^R^-kl ( [tf-T^ [ [ ^-^(l^y,Rdxd,t] dy . 



) 

Denote 



™ = ^i^ JJ -^^y^^dydt. 

Observe that since the solution u is continuous, F : B{0, Rq) M is continuous as 
well. 

Cover B{0,Ro) in n cubic cells, {Ci} of linear size R/2. Note that 

4 < n < 8 

and the area of a cell Ci is 

A(a) = ^ . 

If a cell intersects the sphere iS'(0,i?o), we extend F to the whole cell by setting 
^(y) = on Ci \ B{0,Ro). Naturally, this extension makes F is measurable (but 
not necessarily continuous) on UQ. 

Let e > 0. Since F is bounded, there exist yj,y. G Q such that 

^^(yz)>supF-^ and F(y)<infF+^. 

Consequently, 

^ « ^ n 

e 



B{0,Ro) UCi 



1 f rV^ , , 1 /i?'' 

Note that since F > and F = outside i3(0, i?o), without loss of generality 
we may assume e i?(0, i?o)- Moreover, the balls {i3(yi,i?)} form an optimal 
covering of -6(0, _Ro) in the sense of Definition 13. II with K2 = S^. Thus, 

R^F{y{) ^Y.^ jj ^'/'y^K d^dt <K2^ [[ ^'/'o dxd< = K2Rleo , 

and so 

1 / i? " ' 



B(0.i?o) 

for any e > 0, which implies the upper bound in the first relation in p.22p . 
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To obtain the lower bound, proceed similarly, 

B{0,Rf,) UCi 



>-f-VEny)--f-" 



Note that even if y. ^ B{0,Ro), we still can choose V'y .-R satisfying (|2.9p - (|2.1ip 
and so the supports of V'y .fl will still cover -8(0, Rq) and 



■p. , , 1 /"/"luP 

T jj 2 ~ T yy 2 

Consequently, 

-8(0, flo) 

and, since e > is arbitrary, we obtain the lower bound in the first relation of 
([8:221) . □ 

The lemma above allows us to to note that the estimates for the optimal ensemble 
averages, {■) r — ^ X^iLi ' that will follow will also be valid for the uniform averages, 

4. Enstrophy cascade 

Let {-B(xi, -R)}i=i^„ be an optimal covering of B{0,Ro). 
Note that the local enstrophy equation (|2.5I) and the definitions of Pr and 
( see (ISH) and ([5TT|) ) imply 

(4.1) vl/^ = ^P^_i^li_ ffl\u\^{dtcl), + „A<j>,)dKdt 

^ i=l 

where (pi — rj-ipi and ipi = ip:xL- r. is the spatial cut-off on B{xi,2R) satisfying 
If 

(4.2) T>^, 

V 

then for any < R < R^, 

(4.3) ^ ^ 



hence, 



■^R > yPR - Er. 
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Using p.20p we obtain 
(4.4) 

leading to the following proposition 
Proposition 4.1. 



1 C0K2 



(4.5) 



*i?>Cil^P0 ( 1-C2-^ 



with ci — 1/Ki and C2 = C0K1K2 (provided conditions i3.18^[37T9\) are satisfied). 
Suppose that 



(4.6) 



7 



1/2 

for some < 7 < 1. Then, for any R, (cj /j) tq < R < Rq 



(4.7) 
where 

(4.8) 



> ci(l - "f^)iyEo = co,yiyEo 



Co,' 



ci(l-7') 



1-7^ 



To obtain an upper bound on the averaged modified flux, note that from (|3.20|) . 
Pji < K2P0, and hence, (|4.f p implies 



< i^Pr + ^Er < vK2Pa + vCoK2^ Eq. 

If the condition (|4.6p holds for some < 7 < 1, then it follows that for any R, 
{c2'/Vi)to<R<Ro, 



(4.9) 
where 
(4.10) 



1 



Col' 



C2 



C2 



= K2 



Cl,7 = K2 

Thus we have proved the following. 
Theorem 4.1. Assume that for some < 7 < 1 
(4.11) ao<nRo, 
where 
(4.12) 

Then, for all R, 

1 



1 + 



7 



K1K2 



1 

VC0K1K2 



— ao<R<Ro, 



(4.13) 

the averaged enstrophy flux ^ r satisfies 

(4.14) co,^iyPo ci,yi^Po 
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where 

(4.15) co,^ = i— ^ , ci,^ = if2 [l + -rP^l , 

Ai [ K1K2 

and the average {■) r is computed over a time interval [0,T] with T > Rq/v and 
determined by an optimal covering of B{0,Rq) (i.e., a covering satisfying 13.18\) 
and BJ^} ). 

Remark 4.1. The theorem provides a sufficient condition for the enstrophy cas- 
cade. If (|4.1ip is satisfied, then the averaged enstrophy flux at scales R, through- 
out the inertial range defined by (|4.13p . is oriented inwards (i.e. towards the lower 
scales) and is comparable to the average enstrophy dissipation rate in B{0,Ro). 
Note that the averages are taken over the finite-time intervals [0, T] with T > Rq/i^ 
(see (|4.2p ). This lower bound on the length of the time interval T is consistent 
with the picture of decaying turbulence; namely, small v corresponds to the well- 
developed turbulence which then persists for a longer time and it makes sense to 
average over longer time-intervals. 

Remark 4.2. In the language of turbulence, the condition ()4.1ip simply reads that 
the Kraichnan micro scale computed over the domain in view is smaller than the 
integral scale (diameter of the domain). 
On the other hand, (|4.11[) is equivalent to 

1 /"/" I, .12.25-1 . 7, ^ 7^ n2 1 



-JJ |c.|^0^«-Mxrfi< ^^-^i?^- J J |V®u;|^0odxrfi 

which can be read as a requirement that the time average of a Poincare-like in- 
equality on B{0, 2i?o) is not saturating; this will hold for a variety of flows in the 
regions of intense fluid activity (large gradients). 

Remark 4.3. If wc do not impose the additional assumptions (|2.10p and (12. lip 

for the test functions on the balls B{xi,R) (t B{0,Ro), then the lower bounds for 
^ R in (|4.5p and (|4.14p will hold with P replaced by the time-space average of the 
non-localized in space palinstrophy on B{0,Ro), 

2T 



P' = ^J^ J \^®^\^vd^dt. 



-B(xo,-Ro) 

This is the case because the estimate Pr > P/Ki gets replaced with 

Pr > —P' ■ 

' - Ki 

Remark 4.4. If we integrate the relation (|2.5p over B{0,Rq) (instead of summing 
over the optimal covering) and use Lemma [3.21 the in (|4.14p can be replaced 
with the uniform averaged enstrophy flux at scales R, 



S(0,_Ro) 

with Ki = 2^ and ^^2=4^. 
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Remark 4.5. Proceeding similarly as above, but using the energy balance equation 
()2.3p we can derive a sufficient condition for the forward energy cascade; if for some 
< 7 < 1 we have 

(4.16) To<C7i?o, 
then for all R, 

(4.17) —To<R<Ro, 

C7 

the averaged energy flux $fl satisfies 

(4.18) co^^iyE'o <<!>R< ci^-yvE'^ , 

where the constants are the same as in Theorem 14. II 
Note that for a v in H^iB{0, Rq))^ 

I V «) = / IV «) vP dx = - / V • Av dx < Ivl I AvI 



and thus 



> 



|V(8)v| - |Av| 

If we extend the analogy with Poincare inequalities used in Remark l4.2l to this case, 
then the last relation suggests that the Taylor's length scale tq should dominate 
the Kraichnan's scale ctq for a variety of flows characterized by large gradients, and 
so the sufficient condition for forward energy cascade, ()4.16|) . is potentially more 
restrictive then (|4.11|) . which is consistent with the arguments that in 2D flows 
the inertial range for (forward) energy cascade, if exists, should be much narrower 
then the enstrophy inertial range. This fact was in fact established in the Fourier 
settings in [7]. 

5. Existence of inverse energy cascades 

Assume u is a solution of the NSE (|2.1[) which satisfies no-slip boundary condi- 
tions in some bounded region 57 C M^: 

(5.1) u{t,x)\g^^O foralH>0. 

For simplicity, we consider fl = B(0,D) (although more general domains would be 
acceptable). 
Define 

(5.2) e=i J J ^-^rjd^dt 

[0,2T]xa 

and 

(5.3) E^^ J J iVuprydxdt, 

[0,2T]xO 

the time-averaged energy and enstrophy in 17 (localized in time), and 

(5.4) 1 

the Taylor's length-scale for J7 (here 77 is a function of time satisfying (|2.8p ). 
We assume that there exists 7 > and a length-scale < Rq < D/2 such that 
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Figure 2. Regions of support (V'xo,fl)- 



(5.5) e < 7^i?Qi? or equivalently, r < 7i?o • 

In order to define localized fluxes toward larger scales we introduce the following 
cut-off functions. 

Let 1/2 < (5 < 1. Define 

(5.6) D{x,R) = n\B{ii,R) . 

For an Xq in fl and Rq < R < D/2 define the refined cut-off functions (j) = 
0xo,t,_r(^, x) = ri{t)^lj{x), where rj = r]T{t) is defined in (j2.8p and ^ = V'xo,-r(x) is a 
C°° function on ft which satisfies 

0<4><1, on D{xQ,R), 4> ^0 on B{xo,R- Ro), 

^'■'^ with E^<^ and #L<^. 

Figure [2] illustrates the definition of ip in the case B{xq,R) is entirely contained 
in ri. 

Define the localized energy and enstrophy associated to the outer region D{xo, R) 

(5-8) e^oM-^ J J ^-^^I'J dxdt 

and 

(5.9) E^o,R = ^ jj l^ul'-^xo.fl dxdt , 
as well as the total energy flux 

(5.10) l-xo^fl = ^ // + u • V0xo,fl rfxdt. 

Note that since ■0 can be constructed such that = rfU-ip is oriented along the 
radial directions outside the ball B{xo,R), ^xo.fi can be viewed as the flux out of 
B{yio,R) (i.e. into D{xo,R)) through the layer between the spheres S{xo,R) and 
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5(xo, R — Ro). Additionally, (|2.3|) confirms that Cxq./j tends to increase on average 

in the case $xo,-R > 0. 

To show existence of inverse energy cascade we proceed similarly to section |4l 
Note that (j5.1|) implies that the relation ()2.3|) holds for = 0, and so, rewriting 

it in terms of the quantities defined above yields 

(5.11) l>xo,fl = J^^xo.fl - 77; // ^ {dt4'^oM + '^A^xo.i?) dxdt . 



Using estimates analogous to (|4.3I) we arrive at 



rig 



(5.12) 
provided 

(5.13) T>^. 

V 

If < i?o < ^ < -0/2, we only need two regions I?(xi, R) and £'(x2, i?) to cover 
f2 (by choosing Xi,X2 S with |xi — X2I > 2i?). These regions provide optimal 
covering of in the spirit of Definition 13.11 which will be used in this section. 

For these optimal coverings we have 

(5-14) < Cfl, = i (exi.i?^ + 6x2, fl) < e 

and 

(5.15) \e<Er^\^ {E^,m + E^.,r) < E . 

Thus, if we sum up ()5.1ip over xi and X2 and use (|5.12p . we obtain the following 
bounds on the ensemble average of time-averaged local fluxes at scales R 

(5.16) ^R^l (*xi,i?, + -ix^.i?) <Er + ^en <vE+^ 
and 

Co _ 1 Co 
— n-ep > -vh ^ 

Rq ^ Rq 

Consequently, 



e 



(5.17) > Er - -^CR > -vE - -^e 



(5.18) ""-E - 2Co^ j < ^R <uEiyl + Co^ 

Going back to the (15.51) we obtain the following. 
Theorem 5.1. Assume 

(5.19) T < -yRo 

for some < Rq < D/2 and < 7 < l/-\/2Co- Then, for all R satisfying 

(5.20) R^<R<£ 
we have 

(5.21) cq^^vE < < co,.yiyE 
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where 



(5.22) 



Co, 7 — 




1 + C07' , 



w/iiZe the averages are taken with respect to optimal coverings and over time inter- 
vals T > Rl/v. 

Remark 5.1. The meaning of the theorem above is that if the condition (j5.19p 
is satisfied, then for a range of scales i?, the average backward energy flux is com- 
parable to the total energy dissipation rate vE. Thus we have a backward energy 
cascade over the inertial range defined by (|5.20p . The sufficient condition (|5.19p 
does not call for r to be much smaller then the internal integral scale Rq. However, 
the inertial range for backwards energy cascade is wide provided i?o ^ £'/2, which 
means that backwards energy cascade will exist for a wide range of scales provided 
T <^ D (according to (|5.20p it will start at scales comparable with t and end at 
scales comparable with D). In particular, the scales D and Rq do not have to 
coincide. 

Remark 5.2. By combining Theorem 15.11 with Remark 14.51 we note that if on 
some ball B{xt),Ri) C f2 the local Taylor scale satisfies tq ^ while the global 
Taylor scale t <^ D, we have both inverse energy cascade on over the range of 
scale satisfying (|5.20p as well as the direct energy cascade inside i?(xo, i?i) over the 
range of scale defined by (|4.17p (with Rq replaced by i?i). 

Remark 5.3. Since u is zero on dH., we may replace in Theorem 15. II the ensemble 
average $7? with the uniform space average 



Remark 5.4. We work with the no-slip boundary condition on ft, but the results 
of this section (with slightly modified e and E) will hold for space periodic or 
vanishing at infinity flows as well. 



Let xq e B{0,Rq), < i?2 < i?i < Rq- In order to study the enstrophy 
flux through the shell A{xo, Ri, R2) between the spheres S{xo,R2) and S{xo,Ri) 
we will consider the modified cut-off functions (p = 0xo,T,_Ri,i?2 (^i ^) = vit)ip{'x.) 
to be used in the local enstrophy balance (|2.5p where rj = r^rif) as in (|2.8p and 
ip = V'xo,fli.i?^2 £ 2?(^(xo,2i?i,i?2/2)) satisfying 




6. Locality of the averaged fluxes 



(6.1) 



0<V<V'o, V' = lonA(xo,i?i,i?2)nB(0,i?o) 
IWI ^ Co lAi^l ^ Co 



where ipo is defined in (|2.7p and 



(6.2) 



R — i?(i?i, R2) 



min{i?2, i?i — R2} ■ 
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Use (j) to define the time-averaged energy, enstrophy, and palinstrophy in the 
shell between the spheres S'(xo,i?2) and S'(xo,i?i) by 



TJJ2 



Then, 



(6.4) 



'^Xo,-Rl,i?.2 



0'xo,fli,/?,2 



\ 1/2 

exo,i?i,«2 1 



Pxo,Ri,R2 



are the ZocaZ Taylor and Kraichnan length scales associated with the shell A{kq, i?2)- 

Also define the localized time-averaged flux through the shell between the spheres 
S{xo,R2) and 5(xo,-Ri) as 

(6.5) *xoA,i?,2 = ^ yyil^pu-V^dxdt . 

Note that (j> can be chosen radially (almost radially in case A{ko, R2) ^ -6(0, Rq)) 
so that 5'xn,-Ri,i?2 — *xo,-Ri ~ ^xo,i?.2/2- Moreover, (|2.5I) implies that this flux con- 
tributes to increase Exo.Ri,R2 on. average. Thus ^xo,_Ri,_r,2 can be viewed as total 
enstrophy flux into the shell A(xo, i?2)- 

Similarly, total energy flux into the shell A(xo, -R2) is defined by 

(6.6) <i'xo,fli,fl2 = ^ (^1"!"+^') u-V^dxdi. 

Note that (f> satisfies similar estimates to (|4.3p (with i? replaced by R), and so, 
if T > Rq/i^, the local enstrophy balance ()2.5p leads to 

^'xo,fli,fi2 > vP^a.Ri.R2 - '^^E^o,Ri,R2 

(6.7) / ^2 \ 



i?2 



for any Xq G B{0, Rq) and any < i?2 < -Ri < Ro- 

Similarly, utilizing (j2.5|) again, we obtain an upper bound 



^'xo,fli,_R2 < '^-Pxo,_Ri,i?,2 + ^■E'xo,-Ri,i?2 
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Combining the two bounds on 4'xo,fli,i?2 obtain 
(6.9) 

thus, we have arrived at our first locahty resuh. 

Theorem 6.1. Let < 7 < 1, xq G B(0, Rq) and < R2 < Ri < Ro- If 

(6.10) cr^o,Ri.R2 < —rj2^ 

with R defined by id.S^) . then 

(6.11) (1 - 7^) J^Pxo.fli.fl. < ^-xo.fli.fl. < (1 + 7') '^P^oMiM2 

where the time average is taken over an interval of time [0,T] with T > R^/v- 

Remark 6.1. The theorem states that if the local Kraichnan scale crxo,fli._R2j 
sociated with a shell A(xo, R2), is smaller than the thickness of the shell R (a 
local integral scale) , then the time average of the total enstrophy flux into that shell 
towards its center xq is comparable to the time average of the localized palinstro- 
phy in the shell, Pxo,Ri.R2- Thus, under the assumption ()6.10p the flux through the 
shell A(xo, i?2) depends essentially only on the palinstrophy contained in the 
neighborhood of the shell, regardless of what happens at the other sales, making 
(|6.10p a sufficient condition for the locality of the flux through v4(xo, Ri, i?2)- 

Remark 6.2. Similarly as in the case of condition (|4.11l) . we can observe that 
condition (j6.10p can be viewed as a requirement that the time average of a Poincare- 
like inequality on the shell is not saturating making it plausible in the case of intense 
fluid activity in a neighborhood of the shell. 

In order to further study the locality of the enstrophy flux, we will estimate the 
ensemble averages of the fluxes through the shells A{xi, 2R, R) of thickness R — R. 
Since we are interested in the shells inside i3(0,i?o), we require the lattice points 
Xi to satisfy 

(6.12) B{x„R)cB{0,Ro) . 

To each A{xi,2R, R) we associate a test function <j)i = i-jipi where 77 satisfies ()2.8p 
and ipi satisfies (|6.ip with xq = x^ and R = R. 

If A{xi,2R,R) (/_ B{0,Ro) (i.e. we have B{xi,R) C 5(0, i?o) and B(xi,2i?) \ 
S(0,i?o) 7^ 0), then e V{B{0,2Ro)) with iP, = 1 on A(xo,2i?,i?) n B(0,i?o) 
satisfying, in addition to (|6.ip . the following: 

"01 — 00 on the part of the cone in centered at zero and passing 
through 3(0, Rq) f] B(x,, 2R) between S{0, Rq) and S{0, 2Ro) 



(6.13) 
and 



ipi =0 on B{0, Rq) \ A(xi, 4i?, R/2) and outside the part of the 

(6.14) cone in centered at zero and passing through S'(0, Rq) n i3(x,;, 4i?) 

between S{0,Rq) and S{0,2Rq). 

Figure |3] illustrates the definition of ipi in the case A{xi,2R, R) is not entirely 
contained in B{0,Rq). 
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region containing the rest of supp(^^) 




Figure 3. Regions of supp(?/;i) in the case A(xi, 2i?, R) (/_ i?(0, i?o)- 



Similarly as in the previous section, we consider optimal coverings of B{0,Rq) 
by shells {A{xi,2R, R)}i^i^n such that (|6.12l) is satisfied, 

and 

(6.16) any x € B{0,Ro) is covered by at most K2 shells 74(xi, 4i?, i?/2) . 
Introduce 

1 " 



(6.17) 



62^:,, 



,2R,R 



E- 



2R.R. 



2RM 



- -J2E^,,2R,R i E'2R,R- - ^E^..2R,R ) , 
2=1 \ i=l / 



1 " 

— — / , Pxi.2R,R , 



and 



(6.18) 



1 " 

^2R,R = - V$x,, 

1=1 
1 " 



2J?,K 



2_R,fl 



Xi,2_R,fi 



i=l 



the ensemble averages of the time-averaged energy, enstrophy, palinstrophy, and en- 
ergy and enstrophy fluxes on the shells of thickness R corresponding to the covering 
{A(x„2i?,i?)},=i,„. 
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Taking the ensemble averages in (|2.5|) and applying the bounds for derivatives 
of (t>i, we arrive at 



(6.19) ^2R,R > vP2R,R - E2R,R 

provided T > Rl/v. 



0/ „^-™ 

If the covering is optimal, i.e., if (|6.12[) and (|6. 15116.1^ hold, then 

1 1 / i? ^ ^ 



(6.20) P2^^>_P> P„ 

n Ki\RoJ 

and 

(6.21) E2R,R<—E<K2(^\ Eo 
where 

(6.22) Po = ^ JJ \\/ ^oj\^(j>odxdt = RlPo 

is the time average of the localized palinstrophy on B{0, Rq) and 

(6.23) E, = ^^ll |c.|2^f-irfxd< = RlEo 

is the time average of the localized enstrophy on B{0,Ro) with (po is defined by 
Let us note that 

utilizing (|Og)) . 1^:^ and (15:^ in the inequality ((gl^ gives 
(6.25) 

Taking the ensemble averages in the localized enstrophy equation (|2.5p again, 
this time looking for an upper bound, yields 

Co 

^2R,R < vP2R,R + 1^^E2R,R ■ 

If the covering {A{xi,2R, R)}i^i n of B{0,Ro) is optimal, then, in addition to 

6.26 P2R,R <—Po<K2[ — ] Po; 

n \Ro/ 

hence. 
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Collecting all the bounds on ^2R,r we obtain 
(6.28) 

which readily implies the following theorem. 

Theorem 6.2. Assume that the condition (^TTTp holds for some < 7 < 1. Then, 
for any R satisfying ^.13^ , the ensemble average of the time-averaged enstrophy 
flux into the shells of thickness R, 5'2i?,,_R, satisfies 

(6.29) co,^ (^-^^ lyPo < ^2B..r < ci,^ (^-^^ yPo 

where c, co.-y, and ci^^ are defined in 14-12^ and j-j-lS^ and the average is computed 
over a time interval [0,T] with T > Rq/v and determined by an optimal covering 
{A(x„2i?,i?)}j=i,„ of B{0,Ra) (i.e. satisfying [6l^) . |gJ5)] . and ![6JB\) ). 



Note that if 

denotes the ensemble average of the time- space averaged modified energy flux 
through the shells of thickness R then, dividing (|6.29p by R^ , we obtain the follow- 
ing. 

Corollary 6.1. Under the conditions of the previous theorem, 
(6.30) co^^iyPo < '<i^2RM < cu^iyPo . 



Theorem 16.21 allows us to show locality of the time-averaged modified enstro- 
phy flux under the assumption (|4.1ip . Indeed, the ensemble average of the time- 
averaged flux through the spheres of radius R satisfying (|4.13p is 

= R^^R . 

According to Theorem 14. 1[ 

Co,7 (^-^"j < < ci^-y (^-^^ lypQ . 

On the other hand, the ensemble average of the flux through the shells between 
spheres of radii R2 and 2i?2 , according to Theorem 16.21 is 

co,7 (J^^ - ^^R2,R2 ^ ci,^ ' 

Consequently, 

(6 31) ^^2i2 ( < ^2^2,-B.2 ^ £Vy ( ^ ^ 

ci^^ \ R J ~ §^ ~ Co, 7 V R 

Thus, under the assumption (j4.11l) . throughout the inertial range given by (j4.13l) . 
the contribution of the shells at scales comparable to R is comparable to the total 
flux at scales R, the contribution of the the shells at scales R2 much smaller than R 
becomes negligible (ultraviolet locality) and the flux through the shells at scales R2 
much bigger than R becomes substantially bigger and thus essentially uncorrelated 
to the flux at scales R (infrared locality). 
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Moreover, if we choose R2 = 2^R with k an integer, the relation (|6.31|) becomes 

(6.32) ^2^^ < "^^"t.^^-^"^ < ^22fc , 

Cl,7 Vpfl Co, 7 

which imphes that the aforementioned manifestations of locahty propagate expo- 
nentially in the shell number k. 

In contrast to (|6.3ip . since Eq = RqEq, Pq = RqPq, ^2^2,^2 = ^2 ^2^2,^2 ^^'^ 

= R^^R, 

(6.33) ^ < < ^ , 

Cl,7 Co,7 

i.e., the ensemble averages of the time- space averaged modified fluxes of the flows 
satisfying (|4.1ip are comparable throughout the scales involved in the inertial range 
()4.13p which is consistent with the existence of the enstrophy cascade. 

We conclude this section by noticing that the remarks similar to those at the 
end of section m can be applied here. Namely we have the following. 

Remark 6.3. If the additional assumptions (I6.13|) and (|6.14l) for the test functions 
on the shells A{xi,2R, R) which are not contained entirely in B{0,Ro) are not 
imposed, then the lower bounds in (|6.25p and (I6.29P hold with Pq replaced by the 
time average of the non-localized in space enstrophy on B{0, Rq), 

2T 

Pq^^J J \V ®uj\^vd^dt^ RlP' . 

B(xo,-Ro) 

This is the case because the estimate (I6.20p gets replaced with 

1 f R^^ 



Also, the estimates (|6.31l) and (j6.33p will contain the terms Pq/Po(= Pq/Pq) in 
the lower and Pq/Pq in the upper bounds. 

Remark 6.4. If we integrate the relation (12. 5p over B{0, Rq) (instead of summing 
over the optimal covering) and use Lemma 13.21 the 5'2_R,i?, in Theorem 16.21 can be 
replaced with the uniform averaged enstrophy flux into shells of thickness R, 

1 

•-0 



S(0,flo) 

with Ki = 2^ and isTa =4^. 

Remark 6.5. Working with (|2.3I) yields similar results for the locality of the energy 
fluxes ^xo,Ri,R2 and ^2R,r- Namely, Theorems 16.11 and 16.21 hold with ^' replaced 
with <&, P replaced with E' and length scales a in the sufficient conditions (|6.10p 
and (14.111) replaced with r. 

The locality of energy flux into shells related to the inverse energy cascades is 
established in similar way (except, because of the no-slip boundary condition (jS.ip 
we can set tpo = 1). Note that the flux on a shell is deflned exactly in the same 
way as in (|6.6p . and we obtain the exact equivalent of Theorem 16. II in this setting. 
If in addition the sufficient condition for inverse energy cascade (I5.19P holds, then 
for Rq < R2 < Ri < D /2 we can prove the following equivalent of Theorem [ 
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Theorem 6.3. Assume that the condition i5.19]} holds for some < 7 < l/^2Co. 
Then, for any R satisfying i5.2U\) . the ensemble average of the time-averaged total 

energy flux out of the shells of thickness R, $2i?„i?,j satisfies 

(6.34) "^"'^(^l;) ^ ^2i?,fl < ci,^ (^-|-^ lyE, 

where E is as in h5.S\l . Co.-y and ci_^ are defined in i5.22\) . and the average is 
computed over a time interval [0,T] with T > R^/v and determined by an optimal 
covering {A{x.i,2R, R)}i-i^n of B{0, Rq) (i.e. satisfying i6.1S\) . h6.15]) . and 116. lS\) ). 
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